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We consider the voter model with M states initially in the system. Using generating functions, 
we pose the spectral problem for the Markov transition matrix and solve for all eigenvalues and 
eigenvectors exactly. With this solution, we can find all future probability probability distributions, 
the expected time for the system to condense from M states to M — 1 states, the moments of 
consensus time, the expected local times, and the expected number of states over time. Furthermore, 
when the initial distribution is uniform, such as when M = N, we can find simplified expressions for 
these quantities. In particular, we show that the mean and variance of consensus time for M — N 
is — 1)^ and i(7r^ — 'A){N — 1)^ respectively. 


I. INTRODUCTION 

The voter model is a well studied model in social opin¬ 
ion dynamics [H-Q. In the binary case, each node in a 
network is endowed with one of two states. In a single 
update, a node is chosen randomly and adopts the state 
of a randomly chosen neighbor. Although the model typ¬ 
ically specifies that the states are binary, we study the 
case in which there are initially M states il. We as¬ 
sume M can take any value from 2 to N, where N are the 
number of nodes in the network. As the system evolves, 
it inevitably eliminates a state completely. That is, there 
will almost surely be M — I distinct states in the network 
at some finite future time. This process is repeated un¬ 
til consensus is reached and one opinion dominates the 
network. 

The microscopic rules are unchanged by introducing 
more states in the multi-state extension that we consider. 
Other extensions that modify the update rules by intro¬ 
ducing more states have also been studied The two 
word Naming Game, for instance, is another social and 
linguistic model that updates with different rules than 
the voter model and has noticeably distinct features [ 13 - 
[T3| . Multi-state extensions allow for the consideration of 
neutral opinions, and the behavior of these individuals 
has a vital role in the resulting features of the system. 
For the Naming Game, this intermediate state assumes 
that individuals are flippant in the sense that they al¬ 
ways convey an extreme opinion with equal probability. 
When compared with the three state voter model, the 
resulting features of the model are in great contrast with 
those of the Naming Game [llj, [l3|- The third state in 
the voter model can be thought of as a strongly neutral 
position M- In this context, neutrality itself constitutes 
an opinion of its own that is independent of the other 
states. Moderate individuals who accept a compromise 
between the two extreme beliefs while rejecting the ex¬ 
tremes themselves would fall into this category. They 
would tend to speak the moderate opinion instead of ei¬ 
ther extremes. This paper will provide much insight into 
the properties of social systems with strong neutrals in 
contrast to the weak neutral cases that have been con¬ 


sidered previously in the Naming Game. 

One may observe that many social discussions pose 
a dichotomy between opposing viewpoints. Despite the 
unique characteristics of each individual person, like- 
minded groups often form quickly and dominate the dis¬ 
course. This is certainly evident when considering politi¬ 
cal parties, which show that only a few distinct opinions 
are expressed on a national scale. It may seem counterin¬ 
tuitive that these associations form so quickly when indi¬ 
vidual thought would initially provide for much disagree¬ 
ment. We will show that the multi-state voter model rec¬ 
onciles how social and political systems quickly condense 
into only a few dominating opinions for any value of M. 
We will show that when M = N, the expected time to 
consensus is barely larger than when M = 2. Further¬ 
more, it will be shown that 0{N) opinion states will be 
eliminated in 0(1) time. These results indicate that the 
opinions of individuals quickly condense together into a 
few dominating groups and that the system reaches con¬ 
sensus at a much slower rate. 

We will outline a very powerful procedure that can di¬ 
agonalize the transition matrix for the multi-state voter 
model. This procedure is a generalized application of the 
methods used to diagonalize the binary state voter model 
[l5l |. The model is viewed as an urn problem in which 
two balls are chosen randomly and placed in the urn from 
which the second ball was chosen. By using generating 
functions, we will solve for all eigenvalues and eigenvec¬ 
tors of the transition matrix, which in turn allows us 
to easily find all future probability distributions. These 
techniques can be traced to the solution of the Ehren- 
fest urn model [l3|- The multi-state voter model can be 
thought of as an M urn system will N balls distributed 
amongst them with the same update rules. To solve this 
urn system by generating functions, we cast the spectral 
problem as a partial differential equation that can be 
solved. This the most natural approach since each inde¬ 
pendent variable in the generating function corresponds 
to an urn/opinion state. 

The paper will be outlined as follows. In Sec. im we 
clearly define the complete graph model and the nota¬ 
tion that will be used. In Sec. imi we will consider the 
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Markov transition matrix for a single update and solve 
for its eigenvalues and eigenvectors. We use generating 
functions to solve the spectral problem [l^, [l^. This 
solution has several applications and consequences that 
will be described in detail in Sec. El In particular, we 
consider the m step propagator, moments of consensus 
time, expected local time, expected time to collapse, and 
exact moments of consensus time. Also in Sec. E we 
consider the case when the initial distribution is uniform, 
which is always the case when M = N. In uniform cases, 
the solutions we provide simplify considerably. 


II. THE MULTI-STATE VOTER MODEL 

We assume throughout that there are N nodes and 
that every node is connected to all other nodes. We ap¬ 
ply the model on the complete graph primarily for ana¬ 
lytical tractability. Each node in the network is initially 
endowed with one of M possible opinion states denoted 
by Ai, A 2 ,..., Am- We define the components of n(m) 
to be total number of nodes with opinion Aj. The voter 
model prescribes the random walk for the macro-state 
vector n. That is, we can write for time step m. 


< ni(m -1- 1) ^ 

- 

( ni{m) \ 


( Ani(m) \ 

\nM{m + 1)/ 


\nM{m)J 


\AnM{m)J 


Here, An contains the random nature of the walk at time 
step m. In a single update, Arii = 1 and Arij = — I 
for some i and j, which implies that An = — Sj for 

standard basis vectors . The probability that An takes 
this value is prescribed by the rules of the voter model 
and is given by 


Pr{An = e,-ej|n(TO) = a} = -^^^^^^. (2) 

This accounts for all lazy steps in the system as well since 
there is a non-zero probability that An = 0. 

We define the macro-state probability distribution by 
= Pr{n(m) = a}. We now define a generating 
function for the probability distribution of macro-states 
at time step m as 

QW(=r)= ^ (3) 

|a|=iV 


The vector power is interpreted in the sense of the multi¬ 
index notation of Laurent Schwartz [T^. which shall be 
used extensively. This generating function allows us to 
very easily find a succinct expression for the Markov tran¬ 
sition matrix for a single step of the multi voter model. 
The form of the generating function allows us to de¬ 
termine the shift and differentiation properties of Q^'^\ 
These properties are given by 


1 . 


OiiCtj (m) ^ XiXj /^{‘^) 

N{N-1)^^ ^ 7V(7V-1) 


2 . a 


(m) 

ot — ei-\-ej 




Using these properties, we can rewrite the spectral prob¬ 
lem as an equivalent partial differential equation for the 
generating function of the macro-state probability distri¬ 
bution as 


Qi^+l) 


M-1 M 

0 '"' = E E 

2=1 J = 2+l 


{xi — XjY 
N{N — \) dxidxj ’ 


( 4 ) 


This constitutes a transition matrix that we wish to di¬ 
agonalize. Given the diagonalization of the transition 
matrix, we can find all future macro-state probability 
distributions explicitly, which yield several exact solu¬ 
tions. To accomplish this, we proceed to solve for all of 
its eigenvalues and eigenvectors. 


III. SPECTRAL SOLUTION 

We can solve the partial differential equation given in 
Eqn. (|4]) exactly for all eigenvalues and eigenfunctions. 
To do this, we write the eigenvalue problem in generating 
function form. Eor eigenvalue A with eigenvector v with 
components c^. Let 

G(x) = ^ c„a:" (5) 

|a|=7V 


be the generating function for the eigenvector v. Eur- 
thermore, we require that each component in the vector 
a is non-negative. This is because the index o: directly 
represents the number of individuals with each opinion 
type. We can rewrite the eigenvalue problem for Eqn. 
0 as 


N{N -1){X-1)G = 


M-l M 

i=l j=i+l 


d^G 

dxidxj 


( 6 ) 


We solve for both A and G by utilizing a linear change 
of variables x ^ u and G(x) = H{u). Since the change 
of variables is linear, we expect H to have the same form 
as G. So, we define so that 


Hiu) = 6„u“. (7) 

|a|=7V 


The change of variables is chosen so that the resulting 
difference equation for is explicit. Having an explicit 
equation for ba will allow us to to find all eigenvalues 
and eigenvectors exactly. The change of variables that 
accomplishes this is given to be 
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Ui= Xi- XM (8) 

: (9) 

y-M-l = Xm-1 — Xm (10) 

Um = Xm- (11) 

With this change of variables we can write Eqn. (I6|) as 


N{N-l){\ 


M-l 


1)H= E 

2=1 



M-l 

i=i 


M-l 

j=i+l 


( 12 ) 


To simplify this equation, we use the following identity: 


M-l M-l 

'y ' y ' ('*^2 '^j)HuiUj = 

2=1 j—i-\-l 

M-l . M-l ^ 

^ ^ j . (13) 

i^l ^ j^l ^ 

The proof of this identity is given in Appendix Ap¬ 

plying the identity and canceling like terms reduces Eqn. 
(ff^ to 


M-l 

iV(iV_i)(A-l)i?= ^ 

2=1 




M-l 

y ( ‘^'^i^jHuiUj T ^i^UiUM 

j=i+l 


(14) 


We now rewrite this as a difference equation for the co¬ 
efficients of H. By Eqn. we obtain 


M-l 

N{N-l){X-l)b^ = Y, 
2=1 


ai{ai - l)ba 


M-l 

- Y 2aiaj6„ + (ai-l)(aM + l)fea-e i+eM 

j=2+l 


(15) 


Let w(a) = X)fci ^ “j- We use this to reduce Eqn. (fT^ 
to an explicit form given by 

““ fV(iV-l)(A-1) +u;(a)[«;(«)-!]■ ^ ’ 

Recall that each > 0. 

Observe that if Eqn. (1161) is non-singular for every a, 
then every &<-, = 0. Since this corresponds to the triv¬ 
ial solution to the eigenvalue problem, we discard these 


solutions. Requiring a singularity in Eqn. (1161) implies 
that the eigenvalues for the transition matrix of the M 
state voter model are 




(a) = 1 - 


w(a)['u;(Q:) — 1] 

nIn - 1) 


(17) 


Since many values of a. will yield the same value for 'w{ot), 
it is clear that there will be many repeated eigenvalues. 
Given that w{ol) ranges from 0 to N, the set of eigen¬ 
values is the same for all values of M. However, the 
multiplicities of each eigenvalue will vary with M. 

The components of the eigenvectors can be found by 
transforming back from H{u) to G{x). This will yield a 
relationship between Ca and ba. Using generating func¬ 
tion techniques, this relationship is given by 




l/3|=JV 



(18) 


Here, /3 is a multi-index that has M non-negative com¬ 
ponents, similar to a. The mathematical derivation of 
Eqn. (flHl) is given in Appendix [BJ 


IV. APPLICATIONS 

With the solution to the spectral problem available, we 
can exactly calculate several quantities and estimate oth¬ 
ers. Below we dehne and calculate the expected collapse 
times, the moments of consensus time, the expected lo¬ 
cal times, and the expected number of states over time. 
We also consider the special case when M = N, which is 
when all nodes begin with a distinct opinion. The connec¬ 
tion between the spectral problem and these exact solu¬ 
tions is as follows. With all eigenvalues and eigenvectors, 
we can explicitly diagonalize the Markov transition ma¬ 
trix for the macrostates of the system. The probability of 
achieving each macrostate governs each of the following 
quantities, so having an exact m-step propagator allows 
us to exactly calculate these solutions. 


A. Moments of Collapse Times 

The collapse times, Tfc, are the amount of scaled time, 
m/N, until a state is eliminated from the system. That 
is, if there are k states in the system, then the collapse 
time is the time until only k — 1 states are present. 
Once eliminated, a state can never be reintroduced in 
the system. This process is repeated until each individ¬ 
ual adopts a single consensus state. 

Here, we estimate the collapse times for each k. We do 
this by estimating the probability that k have survived 
by at time m. We call this the survival probability and 
denote it by Skijn). Given that the system is martingale 










and that the solution of the spectral problem is known, 
the survival probability for M states can be bounded as 
follows: 
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where da is the initial distribution expressed in the eigen- 
basis. Let 


5fe(m) = 0{XT). 


(19) 


One can verify that the dominant eigenvalues are Xk by 
considering the eigenvector Ca = 1 and assuming that 
the index ct corresponds to macrostates where k states 
have survived. This eigenvector provides a uniform upper 
bound for the survival probability the system. 

To find the moments of the time to collapse, we use 
this to estimate the probability that the system collapses 
at time m. This is equal to the difference between the 
survival probabilities from time step m — 1 to to. That 
is, by Eqn. (US), the probability of collapse at time to 
is 0[X'^~^(1 — Afe)]. Therefore, we can write the pth 
moments of collapse time as 




= Oipl 


m—1 

N- 1 


fc(fc- 1)_ 

The estimate varies for various initial conditions, but is 
asymptotically correct for all TV, k, and p. We will make 
use of this result when determining the asymptotic be¬ 
havior of the moments of consensus time and the ex¬ 
pected number of states over time. 


B. Moments of Consensus Time 

The consensus time, r, is the amount of scaled time 
until every individual adopts a single opinion. All con¬ 
sensus states are absorbing, so once this state has been 
achieved, all dynamics in the system halt. We can use the 
solution to the spectral problem to find all moments of 
the consensus time. Furthermore, we will use estimates 
to find the asymptotic behavior for large N. 

To find the consensus time, we define to be the 
probability that the system reaches consensus at time to. 
This is equal to the probability that the system has only 
one individual that has a different state than all of the 
others and then adopts the majority opinion. Therefore, 
we write 


( 22 ) 

j=i 

Given the solution of the spectral problem, we can rep¬ 
resent the macrostate probability as 


M M 


,(”i) _ 




(23) 


M M 


S(3 = [v/ 3 ]e, + (JV-l)e,-- 

i=l i=l 
3=Xi 


(24) 


With this, the moments of the consensus time are given 
by 




E[r^] = 

m—1 

^ oo M M 

“vITITITITwm.. 


(25) 


+ (iV-l)e, 


.m^ 


m—1 i—1 j — 1 f3 


(20) 

3=Xi 

(26) 

(21) 

1 P^-Sp 

^ N^{i-Xpy+^ ■ 

(27) 


This is an exact solution for the moments of consensus 
time. Note that the quantity sp depends on the initial 
distribution through da- The eigenvectors, vp can be 
determined component-wise by Eqn. (0. 

The formula given in Eqn. (EZD is exact for all N, 
p, and initial conditions. We now extract asymptotic 
information about the moments of consensus time. In 
particular, observe that the consensus time is the sum of 
all collapse times. Therefore, 


(28) 

M 

(29) 



\-i\=p 


E (;) n 


k=2 


Here, 7 is a vector with components 72 ... 7m- The 
multi-index notation is used to denote the multinomial 
coefficients as well. Also, the collapse times are indepen¬ 
dent random variables. So, when taking the expected 
value of r^, we obtain 


E[t^] = 


E ; 

ItI=p ^ ' 


M 




.7fcl 


(30) 
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Using the estimate in Eqn. (EH, this becomes 


E u no 




7, 


M 




M 

=p\N-p ^ 

ItHp fc=2 


[A'(l - Afc)]^^ 
1 


(1-Afc)T'U 


dependence on M affects an 0(1) factor of the moments 
of consensus time that one may not wish to casually ig¬ 
nore. For instance, when M = 2, we have rj(2,p) = 1 
whereas p{M,p) < 3 for large M. This suggests that the 
uniform estimate can be up to three times as high as the 
exact solution as M changes. Furthermore, for fixed M 
(32) and as p —>■ oo, the estimate given by Eq. dSD is exact. 
Therefore, we have that 


= O 

= o 


M r 


p\N<- Y, n 


“I 7fc 


k{k-l) 


M r 


p\Nn-p Y. n 

lo'hp fc=3 


k{k - 1) 



(33) 

(34) 




We take the big-O outside of the product because as 
the system evolves, the macro-state probability distribu¬ 
tion tends to the uniform distribution, which corresponds 
to the dominant eigenvalue in the system. This means 
that the estimate given in Eqn. m without the big- 
O is the exact solution. Furthermore, the moments are 
bounded by the dynamics when M = N, which exam¬ 
ined in Sec. IIVEI The initial condition also may provide 
further dependence on M, however this dependence is 
bounded, which does not affect the validity of the result. 

Let 


M r 

7i{m,p) = X! n 

\-i\=P fc=3 


-| Ik 


fc(fc-l) 


with 


7 = ( 72 , 73 , ■ • ■,7m)- 
We can therefore write 


(35) 


(36) 



P 


E;[rP] = 0[p!iVP2“Pry(M,p)] (37) 

We will now provide some of the fundamental proper¬ 
ties of r]{M,p). Intuitively, the meaning of r]{M,p) is the 
correction made to the estimate by changing M. First, 
we show that r]{M,p) is bounded above: 


M 


v{M,p) < X! n 

0<^<pk—3 


1 7fc 


>(fc- 1) 


< 


n 

k^3 

M 

n 


1 - 

2 

p+i 

k{k-l) 

1 - 

2 


fe(fc-i) 


fc=3 


k{k-l) 
k{k — 1) — 2 


= 3 


M - 1 
M +l 


(38) 

(39) 

(40) 

(41) 


This result shows that the moments of consensus time 
can be estimated uniformly in M by 0{p\NP2~P). The 


FIG. 1. Simulation of the voter model for N = 100 and M = 
50 is plotted with p. For each p, the simulation is averaged 
over 1,000 runs. Since p is bounded, Eqn. dszi predicts a 
linear relationship with p with slope log(lOO) —log(2) « 3.912. 
The best fit line for the data is given, which has slope 3.951. 

For the first and second moments, evaluations of 
irj{M,p) are also given: 

p(M,1) = 2(^1-7^ (43) 

,(«,2) = |(,^-9) + 2(i-T)UjU + 0(M-*). 

(44) 

Eqn. (H51) shows that the expected consensus time is 
always 0{N) regardless of the number of initial states 
M. We use these particular cases will be used explicitly 
in Sec. IIVEI when studying the exact solutions of the 
moments of the consensus. 

For small M, r]{M,p) can be easily calculated exactly. 
For M = 2,3,4, we have 
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FIG. 2. Simulation of the voter model is plotted with log N. 
Data is averaged over 2,000 runs with p = 2, M = 20. Eqn. 
([571) predicts a linear relationship between the second moment 
and logN with a slope of 2. The best fit line for the data is 
given, which has a slope of 2.094. 


Now, Mq,( 0) = 1 if n(0) = a and equals 0 otherwise. 
Therefore, E[Ma{0)] = aa\ which is given by the initial 
condition. Similarly, AMa{m) = 1 if n{m + 1) = a and 
equals 0 otherwise. The probability that AMa(m) = 1 
is . So, the local time for state n = ct is 

OO 

E[M^{oo)] = (50) 

m—0 

CXD 

= E (51) 

m—0 

We use the diagonalization given in Eqn. (IMl) to com¬ 
pute this. We ignore the terms that have eigenvalue 1 
however because these correspond to consensus states. 
We only consider non-absorbing states when consider¬ 
ing local time. Let E[M] take components E[Ma{oo)]. 
Therefore, the local time reduces to 

mi = i: (52) 

V I-A/3 

The components of M that correspond to consensus 
states are meaningless, as it is understood that when the 
system enters a consensus state, the dynamics halt en¬ 
tirely. The other components are exactly equal to the 
expected local time for their respective macrostates. 


r]{2,p) = 1 

(45) 

.3 .3-P 


(46) 

9 fi-P 

?7(4,p) = 5-3 P + —. 

(47) 


Note that as p ^ oo, these solutions for ri{M,p) exactly 
match the upper bound given in Eqn. gu. The depen¬ 
dence on p for each rj{M^p) always takes the form of an 
exponential attraction to the upper bound in Eqn. dm. 


C. Expected Local Times 

The local time is defined as the amount of scaled time, 
m/N, spent at each macrostate n prior to consensus. 
If Mct{m) is the number of times state n = a. has been 
visited by time m, then one can construct a random walk 
model for each Mq,. That is, we write 

Ma{m + 1) = Ma{m) + AMa{m). (48) 

The expected local time, therefore, is E[Mct{oo)] with 
this notation. Taking the expected value of Eqn. (H5I) 
and summing from m = 0 to m = (X), we get 

OO 

E;[M„(oo)] = E[M„(0)] -b E E[AM^{m)]. (49) 

m—0 



FIG. 3. Example of expected local times for N = 30, and 
M = 3. The initial condition is ni(0) = n2(0) = 10. Most 
of the time is spent on the boundary where one of the states 
had been eliminated. Each macrostate on the boundary has 
nearly equal local time. 
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D. Expected States over Time 

Given the solution to the spectral problem and the 
expected collapse times, finding the expected number of 
states over time, s{t), is straightforward. To do this, we 
sum the collapse times, r^, from k = s + l to k = M. 
Using Eqn. (ED, we show that the time for s states to 
exist in the system is 


k — 1 states is given to be exactly 

We now use this to find the exact number of states over 
time. Recall from Sec. IIVDI that the time to achieve 
s states is the sum of collapse times from fc = s + 1 to 
k = M. Therefore, we obtain 


M 

t= Y. o 


= o 


N 


k{k-l)) 



(53) 

(54) 


M 


E 


fc=s+l 


TV- 1 
k{k-l) 




(57) 

(58) 


Here, t is interpreted as the scaled time m/N. Solving for 
S shows that the expected number of states as a function 
of time is 


= (s +1) 

for a constant rate c. This result is in agreement with 
the literature regarding the multi-state voter model @ . 


E. Ordering Dynamics for Uniform Distributions 
and M = N 

While the above solutions hold for all M, N, p, and 
initial condition , the ordering dynamics of the model 
reduce significantly in the special case where the initial 
condition is uniformly distributed. This is because the 
uniform distribution is also an eigenvector for all M. 
The eigenvalue that corresponds to this eigenvector is Afc 
when there are k distinct opinions in the system. There¬ 
fore, the diagonalization reduces considerably, which al¬ 
lows us to find simplified expressions for the above quan¬ 
tities. 

A special case of a uniformly distributed initial con¬ 
dition is when M = N. This is when each individual 
adopts a unique, personal opinion state prior to global 
discussion. In this case, there is only one possible ini¬ 
tial condition, and therefore uniformly distributed. Also, 
notice that in the next iteration of the model, one state 
will have been eliminated with probability 1. In this time 
step, two individuals will have the same state while the 
others possess distinct states. The probability distribu¬ 
tion at this time step is uniform (constant). That is, 
each state is equally likely to have the two individuals 
than any other state during the first time step. 

We can compute exactly for each k for the uni¬ 
form case. Because the probability distribution of the 
macrostates is an eigenvector, the estimates we calcu¬ 
lated above are exact. In particular, the survival prob¬ 
ability is given to be Sk{m) = A™. Thus, making this 
substitution into the derivation given in Sec. IIV A[ we 
find that the expected time to collapse from k states to 


Therefore, the expected number of states is given to be 

When M = 0{N), this shows that 0{N) states will be 
eliminated in 0(1) time. For example, take t = 1 and 
M = N and observe that s(l) ^ N/2. This shows that 
the system retains only half of its initial number of states 
at t = 1, which corresponds to a sweep of nodes in the 
network. For any t = 0(1), we find that only a fraction of 
the initial number of states remain, so 0{N) states were 
eliminated in this time. This shows that these systems 
quickly converge to 0(1) states relative to the consensus 
time. 

Taking s = 1, the resulting value of t is the expected 
time to reach consensus. Doing so shows that the ex¬ 
pected time to consensus is 


E[t] = 


{N - 1){M - 1) 
M 


(60) 


For M = N, this result also shows that the expected 
number of interactions between individuals until consen¬ 
sus is reached is exactly (A^ — 1)^ and that the consensus 
time as close to N. When M = 2, the expected consen¬ 
sus time is at most A^log2 [2l|, which is not much less 
than the M = N case. 

We also expand the methods in Sec. IIV Bl to find all 
moments of consensus time. By combining the observa¬ 
tion in Eqn. (15(11) with Eqn. (1551) . we can find that 


E[tP] = p\{N - l)P2-Pr]{M,p). (61) 

By using Eqn. (PP. note that the expected time to con¬ 
sensus given in Eqn. (1501) agrees with this result. 

To find the second moment of consensus time, we apply 
Eqn. (PD to Eqn. (15TD to show that 


Up] 





9)+ 1 



1 

3M3 


+ 0 



(62) 


















Fig. |4] features this result. We can also use this to find 
the variance of the consensus time. We combine Eqn. 
(l62)l with the the p = 1 case to find that 


Ear(r) 


{N-lf 



9) 


3M3 


O 


M4 J 


(63) 


This shows that the variance of the of the consensus 
time for uniform distributions does not change much 
with M. Furthermore, taking M = N, the first term 
in the expansion makes for a good estimate, with higher 
order terms being 0{N~^). So, for M = N, we have 
Var{T) ~ i(7r^ — 9){N — 1)^. 



FIG. 4. Simulation data for the second moment of consensus 
time, T 2 , with N = 100 over 30,000 runs is plotted as black 
circles. The exact solution given in Eqn. ((621) is given as the 
blue curve. In addition, the upper bound found by applying 
Eqn. (I 4 TI) to Eqn. (I 6 II 1 is plotted in red. The upper bound 
overestimates the data by a factor of 1.107007 at most. 


V. CONCLUSIONS 

The discussions above have shown many solutions to 
the voter model. Some of these solutions, such as the 
expected time to consensus and the expected states over 
time confirm existing results about the ordering dy¬ 
namics found by other techniques. The other solutions 
require the above spectral analysis to address fully. In 
particular, we found all moments of collapse time, all ex¬ 
pected local times, and provided closed form expressions 
for all moments of the consensus time for uniform ini¬ 
tial distributions. We then used this to find the mean 


and variance of the consensus time and showed that the 
variance converges cubicly in M. 

The results and methods found here show that there is 
potential to solve other general opinion problems. Sim¬ 
ilar multi-state models that may be solvable by these 
techniques include the multi-allelic Moran model and the 
diploid Moran model of Genetic Drift [25l - [^ . The K- 
word Naming Game is one such model that has many 
more than 2 states. The Naming Game with iF-words 
has 2^ — 1 words, which makes numerical solution of the 
ODE dynamics difficult The methods and results 
given here are purely analytical and therefore can circum¬ 
vent those numerical challenges. Compared to the M = 3 
case of the voter model, the Naming Game has a much 
smaller consensus time of O(logiV) [l3,113 compared to 
the 0{N) time of the 3-voter model. The variance of the 
consensus time for the voter model is about 6.944 x lO^ 
for N = 500, which is lower than for the Naming Game, 
whose variance is around 1.1 x 10^ [13 with committed 
fraction of 20%. Even though the Naming Game appears 
to give a higher variance, the result also shows that it is 
linear with N, whereas for the 3-voter, we showed that it 
is quadratic. So, for N large enough, the variance of con¬ 
sensus time for the 3-voter model will eventually become 
much larger. 
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Appendix A: Proof of Identity 


Here we prove the identity given in Eqn. 0 that was 
utilized to solve for all eigenvalues and eigenvectors of 
the multi-state voter model. To begin, consider 

M-lM-l M-1 / i M-1 \ 

i=l j=l i=l \j=l j=i+l ) 

m 

For the first double sum on the right hand side, we inter¬ 
change the sums. Because the sums are dependent, we 
obtain 


M-1 M-1 




M-1 M-1 M-1 M-1 

H 


j=i i=j 


i—l j—i-\-l 


(A2) 
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We relabel i j in the first sum on the right hand side 
and separate the i = j term to obtain. 

M-lM-l M-l / M-l \ 

'y ' 'y ' '^iHuiUj = 'y ^ j UiHuim + 'y ^ UjHu.uj j 
i=i j=i i=i \ j=i+i J 


Here, (/3)a is the multi-index Pochhammer symbol, 
which is defined by (/3 i)qi ■ ■ ■ {PM)aM- Because we de¬ 
fined I a I = TV, the only term in the sum that is non¬ 
zero is when f3 = a. Therefore, D°‘G = a!ca, where 
a! = ai!... um - With this observation, we use the defi¬ 
nition of H to obtain 


M-l M-l 


Rearranging terms in this equation shows that 


+ E E ufHu.u,- (A3) G{x) = H{u{x)) 


M-l 


M-l M-l 




j=i+l 


i=l \ j=l 


= E 

\f3\=N 

= E 

l/3|=iV 


"M-l 


{Xi - XmY' 


2 = 1 
'M-l 0i 

HE 

, 2=1 7i=0 


^/5m 




(A4) 


(B2) 

(B3) 

^M ' 

(B4) 


This concludes the proof of the identity of Eqn. m- 


Appendix B: Calculation of Eigenvector Components 


In this section, we utilize generating function tech¬ 
niques to relate ba to Ca- One strategy is to substitute 
u ^ X into the definition of H{u) and combine all terms 
together. By definition, this must equal G{x) and so the 
resulting coefficients must be Cq. For large M, this be¬ 
comes cumbersome, so we propose a more general means 
of finding the relationship for any M using differentiation 
properties of the generating functions. In particular, for 
multi-index derivative operator where |q:| = N, we 
have 


D^Gix) = cp{f3Uxf^-^. (Bl) 

|/3|=JV 


Simplifying the expression on the right side gives 


\I3\=N 0<7</3 

(B5) 

Now we take 11“ of this equation for |q:| = N. We found 
that on the left side, we get alca- Therefore, we obtain 



a!ca 


E 

jPl^N 


M-l 

n 

. i^l 





(B6) 


Therefore, we find that 


Cq; 


55 bp(-i)°^^-^^ 

|/3|=Ar 



(B7) 


which is the desired relationship stated in Eqn. m- 
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